RU-98-40 
February 1, 2008 



00 
Os 
Os 



The Spectrum and Effective Action of SUSY Gluodynamics f\ 

G. R. Farraif] , G. Gabadadze, M. Schwetz 

Department of Physics and Astronomy, Rutgers University 
Piscataway, New Jersey 08855, USA 



Q_i| We study the low-energy spectrum of SUSY gluodynamics using the generating functional for 

Green's functions of composite fields. Two dual formulations of the generating functional approach 
are given. Masses of the bound states are calculated and mixing patterns are discussed. Mass 
splittings of pure gluonic states, in the case when supersymmetry is softly broken, are consistent 
with predictions of conventional Yang-Mills theory. The results can be tested in lattice simulations 
of the SUSY Yang-Mills model. 
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q ; Supersymmetric gluodynamics, the theory of gluons and gluinos, is an extremely useful 

testing ground for various nonperturbative phenomena occuring in conventional QCD. The 
Os 1 Witten index of the SU(N C ) SUSY gluodynamics equals to N c M. Thus, the ground state 

of the model consists of at least N c different vacua parametrized by the imaginary phase 
of a nonzero gluino condensate ffl, 0. The different vacua are related by discrete Z 2 n c 
transformations of gluino fields. Once one of the N c vacua is chosen, the Z 2 n c symmetry 
group spontaneously breaks down to the Z 2 subgroup. 

In analogy with QCD, one expects that in each of those vacua, the spectrum of the 
model consists of colorless bound states of gluinos and gluons. Among those are: pure 
gluonic bound states (glueballs), gluino-gluino mesons and gluon-gluino composites. These 
states fall into the lowest-spin representations of the N = 1 SUSY algebra written in the 
basis of parity eigenstates P, ||. The masses and mixings of these bound states can be given 
within the effective Lagrangian approach. The effective action for N — 1 SYM was proposed 
by Veneziano and Yankielowicz (VY) ||. The VY action || involves fields for gluino-gluino 
and gluino-gluon bound states. However, it does not include dynamical degrees of freedom 
which would correspond to pure gluonic composites (glueballs). 

At this stage we would like to make a digression and comment on the physical meaning 
of the VY effective action. This is not an effective action in the Wilsonian sense. In ref. pj 
the VY action was constructed as a generating functional of one-particle-irreducible (1PI) 
Green's functions J]7| . That means that the VY action, being written in terms of composite 
colorless fields of SYM theory, can be used to calculate various Green's functions of those 
composite variables. Performing those calculations, however, one is not supposed to take into 
account diagrams with composite fields propagating in virtual loops. Loop effects are already 
included in effective vertices and propagators occuring in the action. The simplest kind of 
Green's function one might be interested in is a two point correlator. As we mentioned 
above, the composite operators entering the VY action are the interpolating fields for the 
bound states of N — 1 SYM theory. Thus, a two point correlator (or simply a propagator) 
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of those fields can be used to determine the mass of the corresponding bound state. Hence, 
the effective action, or more exactly the generating functional of 1PI diagrams, which we 
deal with in this paper can readily be used to deduce masses of composite bound states of 
the theory. In what follows the effective action (and effective Lagrangian) discussed will be 
understood in the sense specified above. 

In this report we describe briefly how one can generalize the VY Lagrangian in order to 
incorporate the glueball states into the description [§, |J. 

The classical action of iV = 1 SYM theory is invariant under chiral, scale and supercon- 
formal transformations. Once quantum effects are taken into account, these symmetries are 
broken by the chiral, scale and superconformal anomalies respectively. Composite operators 
that appear in the expressions for the anomalies can be gathered into a composite chiral 
supermultiplet TrW a W a M 

The effective action ofthe model can be a functional of the superfield S 

S = ^(TrW a W a ) Q = A(y) + y/Wf(y) + 6 2 F(y), 

where the VEV is defined for nonzero value of an external (super)source Q 0. (3(g) stands 
for the SYM beta function which is known exactly 0. The lowest component of the S 
superfield, A, is bilinear in gluino fields and has the quantum numbers of the scalar and 
pseudoscalar gluino-gluino bound states. The fermionic component in S is related to the 
gluino-gluon composite and the F component of the chiral superfield includes operators 

corresponding to both the scalar and pseudoscalar glueballs (G 2 and G^G^ respectively). 

Assuming that the effective action of the model can be written in terms of the single 
superfield S, and requiring also that the effective action respects all the global continuous 
symmetries and reproduces the anomalies of the SYM theory, one derives the Veneziano- 
Yankielowicz effective superpotential || 

W V Y(S) = 7 SlnA (1) 

where 7 = — (N c g/16n 2 f3(g)) > 0, ji stands for the dimensionally transmuted scale of the 
model and e ~ 2.71. If one uses this superpotential along with the simplest Kahler potential 
(S + S) 1 ^ 3 , one finds that no glueball operators are present in the Lagrangian. Indeed, all the 
glueball fields enter through the F components of the superfield S. These components have 
no dynamics and are thus integrated out from the action. We would like to argue below that 
one needs to use a bigger representation of supersymmetry in order to accomodate also the 
glueball degrees of freedom in the effective Lagrangian. 

In order to determine how glueballs can be included in the action let us concentrate our 
attention on the expression for the F field. Using the equation of motion for the gluino field 
and for the D component one getsQ 



F = £ + iQ = ^-[Gl + iG^G, v ). (2) 



As we have already mentioned, the F field appears in the VY action without a kinetic term. 
The term bilinear in the F field is proportional to F + F = X 2 + Q 2 . Besides that, there are 
terms linear in the F field in the expression for the effective action, thus, the F field can 
easily be integrated out by means of its algebraic equations of motion ||. 

In order to reveal subtleties of this procedure let us write down the following relation 

Q = ^uapH^, (3) 



^In general, one is not allowed to use the equation of motion if the VEV of the F field is considered. 
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where H^ ua P is a field strength for a three- form potential C ua/ 3, H^ ua/ 3 = d^C ua p — d u C^ a p — 
d a C u ^/3 — dpCvai±- The C^ ua field itself is defined as a color singlet composite operator of 

colored gluon fields A%, C^ a = ^{Aft v A a a - A^d^A" - A%d v A$ + 2f abc A a ^A b I/ A c a ), with 
fabc being structure constants of the corresponding SU(N C ) gauge group. The right-left 
derivative in this expression acts as AdB = A(dB) — (dA)B 

Using these definitions one finds that the expression bilinear in the F field acquires the 
following form 

17+ 7-i sn2 1 tt2 

The second term in this expression is a kinetic term for the three-form potential C^a. As 
before, the £ field can be integrated out, however one should be careful in dealing with the 
C^ va field. 



In ref. [|TT| it was argued that the three-form field C^ va plays an important role in the 
description of the pseudoscalar glueball. The glueball can be coupled to the QCD rf meson 
by means of the C^ va field [|T2|| . In the case of SYM theory the analog of the rj' meson is 



the gluino-gluino bound state which acquires mass due to the anomaly in the U(1)r current 
within the VY approach. Thus, it is natural to attempt to couple the pseudoscalar glueball 
to the pseudoscalar gluino-gluino bound state within the VY action using the three-form 
potential C^ va . 

To elaborate this approach, let us rewrite the SUSY transformations for the components 
of the S superfield in terms of £ and C^ ua (instead of F and F + ) [|H| 

5 C A = v^, 5 ( V = iV2o%d^A + v/2C(S + ^e^dPCT*), 

The set of fields given above forms an irreducible representation of supersymmetry algebra. 
All these fields can be assigned to a supermultiplet introduced in ref. ||10|| . That supermul- 



tiplet is called a constrained three- form supermultiplet [110]. The easiest way to present this 

Eel. 



multiplet is to introduce the following real tensor superfield U 

U = B + i0z-t6z + ^-6 2 A* + --9 2 A + ^-e^Be^Cr^ + 

lo lo 48 

^ + o*d^ + \e 2 e - <j»d»z) + \eW - <9 2 z?) . (4) 

It is a matter of a straightforward calculation to check that the real superfield U satisfies 
the relation| 

D 2 U = -~S. (5) 

4 v ' 



^The quantity Q can also be expressed through the Chern-Simons current as Q = d^K^. Using this 
equation one can deduce the relation between the Chern-Simons current and the three- form potential C va p, 
these two quantities are Hodge dual to each other: = ^e^ af3 C va p. 

■f Despite a seeming similarity, the tensor multiplet U should not be interpreted as a usual vector multiplet. 
The vector field which might be introduced in this approach as a Hodge dual of the three-form potential 
C^ va would give mass terms with the wrong sign in our approach (see section 2), thus, the actual physical 
variable is the three- form potential Cp, va rather than its dual vector field (the Chern-Simons current). 
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Thus, the real tensor multiplet U, defined by the expression fl5|), includes all the components 
of the chiral supermultiplet S. Consequently, using the relation (M the VY action can be 
rewritten in terms of the bigger multiplet U. In addition, the multiplet has also a scalar B 
and fermion z. We will show below that this allows one to include glueball operators in the 
effective action. 

First, let us notice some features of SUSY transformations of the components of the U 
field. The components which are shared by the tensor multiplet U and the chiral multiplet 
S (namely A, ty, £ and C ) transform among themselves, while other fields (B and 
z) are connected by SUSY rotations to the other four components. Furthermore, one can 
define a "gauge" transformation of the U field as the following shift U — > U + Y, where the 
superfield Y satisfies the relation D 2 Y = 0. It is important to notice that by means of this 
" gauge " transformation one can get rid of the B and z fields in the expression for the U 
multiplet. This is the analog of the Wess-Zumino gauge for the tensor multiplet U. Thus, 
any Lagrangian written in terms of the S field only, if reexpressed in terms of the U field, 
is necessarily invariant under the "gauge" transformation defined above. As a result, the B 
and z fields can always be "gauged" away from that Lagrangian. Thus in order to be able to 
retain the B and z fields as dynamical variables one must include terms in the Lagrangian 
which breaks this "gauge" invariance. The simplest term of this type is the quadratic term 
{U 2 )\d- Once such a term is included in the Lagrangian, the "gauge" symmetry becomes 
explicitly broken and the B and z components of the superfield U survive as dynamical 
variables. 

Let us now apply the U field formalism to the VY action. In the case at hand the chiral 
symmetry is spontaneously broken by the gluino condensate. In terms of the chiral superfield 
this corresponds to the existence of a nonzero VEV of the S field (S) = /i 3 . With that in 
mind the appropriate relation between the U field and the chiral multiplet is 

D 2 U = - 1 -(S -(S)). (6) 

The only result of this modification is that the field A in eq. (H) gets replaced by the quantity 
A- (A). 

Now use the relation @ to write the action in terms of the U field. In order to break the 
"gauge" invariance of the VY action we add a term proportional to U 2 to the VY Lagrangian. 
An appropriate term with zero R-charge and correct dimensionality is 



U 2 



(s+sy/ 3 



ID- 



(7) 



Below, we are going to show that once this term is added to the VY action, the following 
fields become dynamical: 

• The B field propagates and it represents one massive real scalar degree of freedom 
(identified later with the scalar glueball). 

• The three-form potential C^ uon which becomes massive, also propagates. It represents 
one physical degree of freedom (identified with the pseudoscalar glueball). 

• The complex field A, being decomposed into parity eigenstates, describes the massive 
gluino-gluino scalar, s, and pseudoscalar, p, mesons. 

• z and \& describe the massive gluino-gluon fermionic bound states. 
Relations between masses of these states will be given in the next section. 
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Based on the arguments given abowe one can write down the effective Lagrangian for the 
lowest-spin multiplets of the N — 1 SUSY YM theory in the following form 



^ = + 7 [(Slog ^-S)\ F + h.c] + \ \v, (8) 

where a and 5 are arbitrary positive constants. One obtains the VY Lagrangian in the 
limit 5 —* oo. In general, higher powers of U can also be added to this Lagrangian. Those 
terms would introduce new quartic, quintic and other higher interaction terms. However, 
the quadratic part of the action which defines two-point Green's functions and masses will 
not be affected. In that respect, the effective Lagrangian (§) can be considered as the one 
describing small perturbations of fields about a vacuum state. 

Let us determine the SUSY vacuum state defined by the Lagrangian The potential 
for the model is a complicated function of the variables present in the U superfield. After 
integration over the auxiliary £ field, the bosonic part of the potential is 



2 



_ 2 |0| 6 + / i 6 -2/i 3 |0| 3 cos3p 3 

"0 r/i Mi i no c-/ An\o i / 1 2 



5(16) 2 |0| 2 5(48 



9a\6\ 4 1 



,3 U|2\ 2 



(B 2u 3 
mW {1 + w cos3p} ~ 3llogl v) ' (9) 

where we introduced the notations <p = an d p = arg</>. 

In order to find the vacuum state one should find the absolute minimum of the potential 
(H). Since we are dealing with a supersymmetric model, the value of the potential in that 
minimum has to be zero. As a result of Lorentz invariance, the VEV of the three-form field 
is zero, i.e. (C Ml/T ) = 0. The VEV of Q is also zero due to the CP invariance of the model. 
After some algebra one finds that the only global, CP invariant minimum of the potential 
(|S[) is given by: (0) = p, (B) = (C) = (p) = 0. The effective Lagrangian (g) describes small 
perturbations of fields about the vacuum state defined by these VEV's. 

We would like to make a comment here. The singularity in the potential at S\4>\ 4 = aB 
indicates that for large perturbations the higher dimensional terms omitted in (§) should 
become important. The same multiplier 1 — f cti appears in front of kinetic terms for the 

scalar fields, so the physical potential is always positively defined. As we mentioned above, 
we are mainly concerned with the mass spectrum of the model which can be studied using 
small perturbations about the ground state, so that the approximation given in (|8|) is good 
enough for our goals. 

The actual physical states described by the action are mixed states. Below, we deduce 
the masses of these mixed physical states. Let us write down the mass and mixing terms of 
the Lagrangian separately. One finds the following pairs of bosonic variables being mixed 
with one another 

__, 9a 2 2 81 2 2 2 2 9a 2n 2 27V2 fa , n 

B - s system : s + -a 7 P s + j^—^p B - ——^- an/ fi Bs; 

C -p system : ^7^2^ P + C ^ + -^-ajppe^d^C ; (10) 



In order to find the physical masses one must diagonalize the corresponding mass ma- 
trices. Concentrate, for instance, on the first row of these expressions which describes the 
mixed state of scalar meson, s, and scalar meson- glueball, B. The former gets mass both 
from the superpotential and £7 2 -term while the latter gets mass only from the f/ 2 -term. 
When the mixing term is switched on, the initially heavier state (s) gets even heavier, and 
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the initially lighter state (B) becomes even lighter than they were originally. Performing the 
diagonalization, one finds that the physical eigenstates are mixed states with the following 
mass eigenvalues 



1 2 9a 2 81 2 2 2 81 2 2 2 L la 1 



2 M - = + T° 7 " ± ~4~ y + 2887 <"> 

Here, the subscript "+" refers to the heavier state s which, without mixing, would have been 
a pure gluino-gluino bound state (the s particle). The subscript "— " refers to the lighter 

state B (B in the absence of mixing). 

Studying the potential of the model, we find that the physical eigenstates fall into the two 
different ''multiplets". Neither of them contain pure gluino-gluino, gluino-gluon or gluon- 
gluon bound states. Instead, the physical excitations are mixed states of these composites. 
The heavier set of states contains: 

• A pseudoscalar meson, which without mixing reduces to the _+ (I = 0) gluino-gluino 
bound state (the analog of the QCD r{ meson). 

• A scalar meson that without mixing is an ++ (1 = 1) gluino-gluino excitation. 

• A mixed fermionic gluino-gluon bound state. 

These heavier states become the chiral supermultiplet described by the VY action in the 
limit that the additional term we have added to the effective Lagrangian is removed. The 
new states which appear as a result of our generalization forms a lighter multiplet: 

• A scalar meson, which for small mixing becomes a ++ (I = 0) glueball. 

• A pseudoscalar state, which for small mixing is identified as a h (I = 1) glueball. 

• A mixed fermionic gluino-gluon bound state. 

We call the reader's attention to an interesting feature of the effective action intro- 
duced here. Although the physical states fall into multiplets whose J p quantum numbers 
correspond to two chiral super multiplets, the action is written in terms of one real tensor 
supermultiplet U. The natural question arises whether the whole action can be rewritten 
in terns of two different chiral multiplets. The relation between a real tensor and chiral 
supermultiplets (the so called chiral-linear duality) was established in ref. |10|. For SYM 



theory the chiral-linear duality was used in refs. |L3|]. Applied to our problem the results of 
refs. 



new c 



10| and [13] can be stated as follows. One introduces into the effective Lagrangian a 
liral superfield, let us denote it by x 

x (y, 6) = cf> x (y) + V26* x (y) + 9 2 F x (y). (12) 

One can find an effective Lagrangian written in terms of two chiral superfields, S and x 
which is equivalent to the expression given in (H). In our case 



c = -(s + sy/ 3 \ D + ~(s+s)^( x + x + f 



a 

S 



+ 

D 



1 



j(S log^ - S) + - X (S - /i 3 ) „ + h.c. . (13) 



16 



F 



Comparing this expression to the VY Lagrangian one notices that both the Kahler potential 
and the superpotential are modified by new terms. The multiplets S and x are independent. 
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We would like to relate this expression to the Lagrangian of the theory written in terms 
of the U field (|8|). If the U field is postulated as a fundamental degree of freedom, then the 

S field is a derivative superfield S 
be rewritten as 



/i 3 — AD 2 U. Using this relation the Lagrangian (p~3|) can 



C 



-{S + sy/ 3 

a 

7 (S log -3 - 



D 



+ r 



D 



S) 



+ h.c. 



+ u(x + x A 



+ 



D 



(14) 



This expression depends on two superfields U and x (S is expressed through U in accordance 
with ([])). However, the dependence on the chiral superfield x is trivial, the combination 
X + X + can be integrated out from the Lagrangian ([14]). As a result one derives 



X + X 



2U 



5(s+sy/ 3 ' 



(15) 



Substituting this expression back into the Lagrangian ( PI) one arrives at the original expres- 
sion (||) where the S field is a derivative field satisfying the relation (Qj. 

Let us stress again that the descriptions in terms of the Lagrangian (|8|) and ( |T3| ) are 
equivalent on the mass-shell. In the Lagrangian (S) the dynamical degrees of freedom are 
assigned to the only superfield U, while in the Lagrangian QT3D the physical degrees of 
freedom are found as components of two chiral supermultiplets o and X- The peculiarity of 
the expression (14) is that the chiral superfield x enters only through the real combination 
X J rX + ■ That is why it was possible to formulate the action in terms only of the real superfield 
U. It is essential from a physical point of view since the component glueball field must be 
real. 

In order to make contact with the results of lattice simulations of SUSY YM model [14 
one needs to consider the model with soft SUSY breaking term introduced via gluino mass. 

The potential of the softly broken model consists of two parts , Vq defined in (|9]) and an 
additional SUSY breaking term H 



V = V + m x Re(^- ( p x + 7 3 ) 



16 



(16) 



where m A = m x . 

One calculates minima of the full scalar potential V. Explicit though tedious calculations 
yield the following results. The VEV of the field does not get shifted when the soft SUSY 
breaking terms are introduced. Thus, even in the broken theory (0) = //. However, the <p x 
(and B ) fields acquire nonzero VEVs in the broken case 



c) 
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9a p, 



and (B) = - 

9a 



(17) 



The shift of the vacuum energy causes the spectrum of the model to be also rearranged. 
Explicit calculations of the masses of all lowest-spin states yield the following results 



M 2 p _ scalar± = Ml - -a7/im A (l±vTT^J(4 ± -j— J, (20) 

where M\ denote the masses in the theory with unbroken SUSY || and x = t^at- In 

these expressions the plus sign refers to the heavier supermultiplet and the minus sign to the 
lighter set of states. One can verify that these values satisfy the mass sum rule to leading 
order in 0(m\): 



£ (-l) 2 ^ (2j + l)M| = 0, 

j 



where the summation goes over the spin j of particles in the supermultiplet. 

Let us discuss the mass shifts given in eqs. ( p~8| - pOD . Consider the light supermultiplet. In 
accordance with eqs. ( JT8| - pOD , the masses in the light multiplet are increased in the broken 
theory. The biggest mass shift is found in the pseudoscalar channel. The smallest shift is 
observed in the scalar channel. The fermion mass falls in between these two meson states. 
Thus, the lightest state in the spectrum of the model is the particle which without mixing 
would have been the scalar glueball. There is a fermion state above that scalar. Finally, the 
pseudoscalar glueball is heavier than those two states. 

Let us now turn to the heavy supermultiplet. In the broken theory the masses in that 
multiplet get pulled down. However, all states of the heavy multiplet are still heavier than 
any state of the light multiplet in the domain of validity of our approximations. The ordering 
of the states in the heavy supermultiplet is just the opposite as in the light supermultiplet: 
the lightest state is the pseudoscalar meson, the heaviest is the scalar, and the fermion, as 
required, falls between them. The qualitative features of the spectrum are shown in fig. 1. 




Figure 1: Qualitative behavior of mass spectrum when passing from SYM to softly broken 
model. 



It is not surprising that the lowest mass state obtained in ( p!8| - |20| ) is a scalar particle. 
This is in agreement with the result of ref. [I5| where it was shown that the mass of the 
lightest state which couples to the operator Gt v is less than the mass of the lightest state 

that couples to GG, in p ure Yang-Mills theory. As a result, the lightest glueball turns out 



to be the scalar glueball One can apply the method of ref. [15| to the SYM theory as 
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well. Due to the positivity of the gluino determinant (see ref. [|T(|) one also deduces that the 
lightest state in softly broken SYM spectrum should be particle. The pseudoscalar 

of that multiplet is therefore heavier. 

Our result that the multiplet containing glueballs is split in such a way that the scalar 
is lighter than the pseudoscalar, and vice versa for the multiplet containing gluino-gluino 
bound states, is consistent with expectations from quark-model lore. In ordinary mesons 
the I = 1 states are heavier than their I = counterparts and the / = gluino-gluino bound 
state is a pseudoscalar, while an / = gluon-gluon bound state is a scalar. It is interesting 
that in SYM with massless gluinos the / = and I = 1 bound states are degenerate, but 
when the gluino masses are turned on one recovers the expected ordering seen in qq states. 

Summarizing, we have shown that the generalized VY effective action can be written 
in two different ways. In one case the fundamental superfield upon which the action is 
constructed is the real tensor superfield U. In another approach all degrees of freedom of 
the model are described by two chiral superfields x an d S. In both cases the spectrum 
consists of two multiplets which are not degenerate in masses even when SUSY is unbroken. 
The spin-parity quantum numbers of these multiplets are identical to those of certain chiral 
supermultiplets. 

We introduced a soft SUSY breaking term in the Lagrangian of the N — 1 SUSY Yang- 
Mills model. The spurion method was used to calculate the corresponding soft SUSY break- 
ing terms in the generalized VY Lagrangian. These soft breaking terms cause a shift of the 
vacuum energy of the model. The physical eigenstates which are degenerate in the SUSY 
limit, are split when SUSY breaking is introduced. We studied these mass splittings in 
detail. We have confirmed that the spectrum of the broken theory is in agreement with 
some low-energy theorems namely the scalar glueball turns out to be lighter than the 
pseudoscalar one. The results of the present pap er can be directly tested in lattice studies 
of N = 1 supersymmetric Yang-Mills theory [14[ . 



Even when SUSY is unbroken, the physical mass eigenstates are not pure gluon-gluon, 
gluon-gluino or gluino-gluino composites; rather, the physical particles are mixtures of them. 
The multiplet which without mixing would have been the glueball multiplet is lighter. As 
a result, those states cannot be decoupled from the effective Lagrangian. This means that 
comparisons of lattice results to analytic predictions based on the original VY action are not 
justified. 

The work was partially supported by grant No. NSF-PHY-94-23002. 
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